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Deconfined regions in relativistic heavy ion collisions are limited to small volumes surrounded 
by a confined exterior. Here the geometry of a double layered torus is discussed, which allows 
for different temperatures in its two layers. This geometry enables one to approach the QCD 
continuum limit for small deconfined volumes with confined exteriors in a more realistic fashion 
than by using periodic boundary conditions. Preliminary data from a study for pure SU(3) lattice 
gauge theory support a substantial increase in a pseudo transition temperature. 
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1. Equilibrium simulations of the deconfinement transition and lab plasma 

Statistical properties of a quantum system with Hamiltonian H in n continuum volume V , 
which is in equilibrium with a heatbath at physical temperature T, are determined by the partition 
function 

Z{T,V) =Tre-"/^ = '£{<l>\e-"/^\<l>), d-D 

where the sum extends over all states and the Boltzmann constant is set to one. Imposing periodic 
boundary conditions (PBC) in Euclidean time T and bounds of integration from to 1 /r, one can 
rewrite the partition function in the path integral representation: 

Z{T,V) = jD(l>cxpl^-J^^'^dzLE{(l>,^)Y (1.2) 

Nothing in this formulation requires to carry out the infinite volume limit. 

Past LGT simulations of the deconfining transition focused primarily on boundary conditions 
(BC), which are favorable for reaching the infinite volume quantum continuum limit (thermody- 
namic limit of the textbooks) quickly with temperature and volume of the system given by 

r = 4r = 7-' N,/N,^oo, N.^oo, Lrfinite, (1.3) 
aN-c Lt: 

where a is the lattice spacing. These are PBC in the spatial volume V = (aNs)^. For the decon- 
finement phase created in a heavy ion colfision the infinite volume limit does not apply. Instead we 
have to take the finite volume continuum limit 

Af,/Ar^ = finite, A^^^oc, finite, (1.4) 

and PBC are incorrect because the outside is in the confined phase at low temperature. E.g., at the 
BNL RHIC one expects to create an ensemble of differently shaped and sized deconfined volumes. 
The largest volumes are those encountered in central collisions. A rough estimate of their size is 

n X (0.6 X Au radius)-^ x c x (expansion time) (1.5) 
= (55 fermi'^) x (a few fermi) 

where c is the speed of light. Here we want to estimate finite volume corrections for pure SU(3) 
and focus on the continuum limit for 

L, = aNs = (5 - 10) fermi . (1.6) 

In the following we set the physical scale by 

r'^ = 174MeV, (1.7) 

which is in the range of QCD estimates with two light flavor quarks, implying for the temporal 
extension 

Li,=aA^^ = 1.13 fermi. (1.8) 
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2. Simulations with cold boundary conditions 

We consider difficulties and effects encountered when one equilibrates a hot volume with cold 
boundaries by means of Monte Carlo (MC) simulations for which the updating process provides 
the equilibrium. We use the single plaquette Wilson action on a 4D hypercubic lattice. Numerical 
evidence shows that SU(3) lattice gauge theory exhibits a weakly first-order deconfining phase tran- 
sition at some coupling lif{N^) = 6/gj{Nt). The scaling behavior of the deconfining temperature 
is 

T' = ctAl (2.1) 

where the lambda lattice scale 

aA, = hm =X{g^) M^'^'^'^^^ .-V(2.og^) ^ (2.2) 

has been determined in the literature. The coefficients Bq and b\ are perturbatively obtained from 
the renormalization group equation, 

11 3 34 / 3 ^ ^ 



Higher perturbative and non-perturbative corrections are parametrized in [|l]] by 

= l+a^e-"'-''^'' + a^g^ + a^g'^ with (2.4) 

ai = 71553750, a2 = 19.48099, as = -0.03772473, a^ = 0.5089052. In the region accessible by 
MC simulations this parametrization is perfectly consistent with an independent earlier one in 
but has the advantage to reduce for ^ to the perturbative limit. 

Imagine an almost infinite space volume V = and a smaller sub- volume V\ = L^. y The 
complement to Vi in V will be called Vq (outside world). The number of temporal lattice links A^^ 
is the same for both volumes. We like to find parameters so that scaling holds in both volumes, 
while Vi is at temperature Ti = and Vq at room temperature Tq. We denote the coupling by jSq 
for plaquettes in Vq and by j8f for plaquettes in Vi . For that purpose any plaquette touching a site 
in Vi is considered to be in V\. We would like to have jSf in the scaling region, say jSf = 6. The 
relation 



^10 ^ 7i ^ «o ^ /A(i8o^)AL ^ fM 



iQi" « ^ = = 'rrM: = -9^, (2-5) 



drives jSq out of the scaling region and (after moving over to strong coupling relations) practically 
to jSq = 0, which we call disorder wall BC. In the disorder wall approximation of the cold exterior 
we can simply omit contributions from plaquettes that involve links through the boundary. MC 
simulations, which we quickly summarize here, were performed in Ref.[^]. Due to the use of the 
strong coupling limit for the outside volume, scaling of the results is not obvious. 
We use the maxima of the Polyakov loop susceptibility 

Zmax = ^[(|/'P)-(|P|)']„,,, P = ZP, (2.6) 
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Figure 2: Estimates of finite volume corrections to the width of the Polyakov loop susceptibility. 



to define pseudo-transition couplings ppf{Ns;N-T:). Our BC introduce an order disturbance, so 
that 

P^pM;Nr)=Pf{N,)+ai^^+ai (^)'+ ■■• (2-7) 

holds. Fits of pseudo-transition coupling constant values lead to estimates of finite volume correc- 
tions to Tc as given in Fig. |T} It is seen from the figure that they are consistent with scaling as the 
curves for =4 and = 6 fall on top of one another. There are no free parameters involved at 



4 



Deconfining phase transition on a double-layered torus 



Bernd A. Berg 





L 


\ 




r 


\ 










1 






2 








\ 

\ 










I 


3 






X 





Figure 3: Double layered torus in two dimensions. 

this point, because the non-perturbative parametrization of the SU(3) lambda scale was determined 
in independent literature. Estimates of finite volume corrections to the full width of Polyakov loop 
susceptibilities at 2/3 maximum are shown in Fig. ^. Again, the curves are found to be consistent 
with scaling. 

Results so far show that for volumes of BNL RHIC size the magnitudes of Tc and width 
corrections are comparable to those obtained by including quarks into pure SU(3) calculations (r^ 
in the opposite, width in the same direction). Similar corrections are expected for the equation of 
state. Previous QCD calculations at finite temperatures and densities should therefore be extended 
to cold BC. 

However, there remain problematic questions about disorder wall BC. Although the results 
show scaling, it is unsatisfactory that the disorder wall BC do not reflect a physical outside volume. 
Two properties are desirable: 

1. Inside and outside volumes are kept in the scaling region. 

2. The spatial lattice spacing a., is the same on both sides of the boundary. 

In such simulations one may again keep Tinside at Tc and study its dependence on Toutside- As outlined 
in the next section, this can be done in the newly introduced geometry [Q] of a double layered torus 
(DLT), but has then to be done for a rather small Toutside interval 

routsidee[7;-Ar,r,], Ar>o. (2.8) 

3. Simulations on a double layered torus 

For the DLT the boundaries are glued together as indicated by the arrows in Fig. |3[ Note 
that interchanging labels 3 and 4 on one of the lattices leads to a situation in which some sites are 
connected by two links and the different geometry of a sphere. 
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With DLT BC in the spacelike directions and PBC in the fourth direction one can simulate at 
two temperatures and each volume becomes the outside world of the other. The two temperatures 
are adjusted by tuning the coupling constants in the two volumes. MPI Fortran code for SU(3) 
simulations on a DLT is given in Ref.[^, ^]. Preliminary numerical results with both temperatures 
in the SU(3) scaling region are compiled in the following. 

Figure Q shows the reweighting in j3i of Polyakov loop susceptibilities on a 8^4 lattice, each 
simulation point corresponding to a different pair of coupling constant values (jSq , jSj"), adjusted so 
that is close to the pseudocritical point. 

Using scaling relations the inverse physical pseudocritical temperatures, l/Tf°^ from all our 
lattice sizes, with Tj'"'"' corresponding to the maxima of the Polyakov loop susceptibilities, are plot- 
ted in Fig. ^ versus the outside temperature Tout. Even for the small range of outside temperatures 
in the SU(3) scaling region, one sees already sizable corrections of T™^". The same data are plotted 
in Fig. ^, extending the Tout down to zero, so that the estimates from Ref . ^ can be included (on 
the 8^4 lattice not transition was found with disorder wall BC as indicated here by l/T™^" = 0). 

In summary, our preliminary DLT results in the scaling region are consistent with the disorder 
wall finite volume Tc estimates. These simulations keep inside and outside temperatures in the 
SU(3) scaling region. To achieve also a continuous spacelike lattice spacing across the boundary, 
one has to use asymmetric couplings in space and time directions. Considerable work remain to be 
done to obtain an overall convincing description. 
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41022, DE-FC02-06ER-41439, by NSF grant 0555397, and by the German Humboldt Foundation. 
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Figure 5: Inverse pseudo transition temperatures versus outside temperatures. 
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Figure 6: Results of the previous figure together with disorder wall estimates. 
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